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By studying varies dynamical processes, including coupled maps, cellular automata and coupled 
differential equations, on five different kinds of known networks, we found a positive relation between 
signal correlation and node's degree. Thus a method of identifying influential nodes in dynamical 
systems is proposed, its validity is studied, and potential applications on real systems are discussed. 
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Dynamical processes on networks have been highly 
concerned 0, 0, EJ > the network topological structure 
often plays a crucial role in determining the system's dy- 
namical features. Various dynamical systems on complex 
networks have been widely investigated, such as transient 
activation in neuron networks || , synchronization 0, Q , 
virus spreading or information diffusion y|. Due to the 
complex connectivity and the dynamical nonlinearity, an 
exact mathematical analysis of general dynamical pro- 
cesses on networks is often highly intricate. In general 
dynamical systems on complex networks, the collective 
behavior of nodes remains unclear. 

Recently, studying dynamical signals (the time series) 
presents a different way to understand the dynamical pro- 
cesses on complex networks. For example, Rcf. 9] shows 
cffcctivity to simultaneously characterize the dynamics of 
several thousands nodes by scaling of fluctuations; Ref. 
[Io| gives an example to reconstruct the nonlinear dy- 
namics in a network from observed dynamical signals. 
Beside theoretical interests, studying the dynamical sig- 
nal may provide practical applications, i.e., for some dy- 
namical systems, it would be difficult or even impossi- 
ble to directly obtain the topological information of the 
background networks, such as financial systems , the 
neuronal functional networks 0], e tc-j thus a key prob- 
lem arises that how to obtain the topological information 
from the dynamical signals. 

In the paper, based on principal component analy- 
sis ^]|, we have investigated the cross-correlation ma- 
trix of the dynamical signals, which reflects the inter- 
dependence of the nodes on the network, and studied 
the relation between the cross-correlation matrix and the 
topological connectivity. We found that in most cases 
there is a positive relation between the nodes' degrees 
and the components of the principal eigenvector, the 
eigenvector corresponded with the largest eigenvalue, of 
the cross-correlation matrix, i.e., the larger degree often 
corresponds to larger components. Thus, we suggest a 
method to obtain the topology from the dynamical sig- 
nals. Three types of coupled dynamical models, namely, 
coupled maps, cellular automata and coupled differential 



equations, and five kinds of complex networks are stud- 
ied. To be specific, the dynamical models are coupled 
logistic maps, avalanche processes [lj, [15| , and integrate- 
and-fire models ^(|. The complex network models in- 
clude the Erdos-Renyi (ER) random graph 0] (which 
has a Poisson degree distribution), the Barabasi- Albert 
(BA) model [13] (which has a power law degree distri- 
bution with exponent A = 3.0), the generalized random 
graphs with scale- free (SF) degree distribution [19|, the 
lattice embedded SF model |2(|, and the nearest neigh- 
bor growth model [ilj (which has an exponential degree 
distribution). In the simulation, the exponent of degree 
distribution for SF random graph and lattice embedded 
SF networks is fixed at A = 3.0 to make the results com- 
parable with that of BA model. 

We start with a general expression of dynamics on a 
network. The connecting network is described by its ad- 
jacent matrix A: Aij equals to I if node i and node j have 
a common edge (then j is said to be a neighbor of i and 
vice versa) and if not, and fcj = X)j=i Aij i s the degree 
of node i, Xi(t) is the output signal of node i. The cross- 



correlation matrix is determined as CV 



(x i x j )-(x i )(x j ) 



where <Ji is the variance: <7j = y (a; 2 ) — {xi) , and the 
average is over a period of length Ld- Here we concern 
only about the strength of the correlation between the 
nodes, rather than the signs, thus we define the corre- 



lation strength matrix: CV 



CVj|. Then its principal 



eigenvector could be determined, and the relation be- 
tween the components of the principal eigenvector and 
the nodes' degrees is investigated in detail. 

Coupled logistic maps. The dynamic of an individual 
node i coupled with its neighbors is described by: 



„*+i 



N 



(1) 



where f(x) is the logistic map f(x) = I — ax 2 , e £ [0, 1] is 
the coupling strength. For a given network, the parame- 
ter space (e, a) has two extreme cases: the synchroniza- 
tion regime and complete non-synchronization regime. In 
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the synchronization region, parts of or all the nodes have 
the consistent motion. The cross-correlation matrix may 
reflect the synchronization motion. Here, we investigate 
the system running in the state which is far from syn- 
chronization. Figure ^ is a typical show of the connec- 
tion matrix and the correlation strength matrix, since the 
states have no particular orderliness, the latter is highly 
noised. The same holds for other network models and for 
other dynamical processes. 




in calculating correlation is Ld — 10 5 , and each group 
of data has been averaged over 300 network ensembles. 
The calculation of correlation begins after 10000 time 
steps, while the transition usually takes a few hundred 
time steps. It should be noted that in the subgraphs the 
ranges of x-axes are not the same due to the different 
degree distributions (SF, exponential, and Poissonian). 
Although there are small deviations and the data scatter 
sparsely for SF random graphs and lattice embedded SF 
networks, the positive relation between the components 
and the degrees is clear. The parameters in our study 
are chosen to satisfy that the dynamics are chaotic, and 
the length of the data for calculating the correlation is 
long enough to eliminate some unpredictable factors. In 
order to achieve this, we have studied coupled maps on 
ER random graphs in detail, compared varies data length 
and coupling strength. As Fig. |3 shows, the parameters 
used in our study is appropriate. Similar studies on other 
network models confirm that the relation is stable under 
variations of dynamical parameters. 



FIG. 1: (Color online) A direct show of network connection 
and correlation strength of coupled maps on a ER random 
graph, with N = 100, (k) — 6. Left: Network connection, 
dot at indicates the connection of node i and j. Right: 

Correlation strength, the brighter, the stronger, a = 1.9, 
£ = 0.9. 
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FIG. 2: (Color online) Components of principal eigenvector 
of correlation strength matrices of coupled maps on networks 
vs. node degree for various network models, (a) SF random 
graph; (b) BA model; (c) lattice embedded SF model; (d) 
Nearest neighbor growth model; (e) ER random graph; (f) A 
comparison of all the network models, same symbols represent 
the same networks in (a)-(e). 

The relationship of the components of the principal 
eigenvector of the correlation strength matrix and the 
nodes' degrees is shown in Fig. [21 In the simulation, the 
degree distribution exponent for SF networks is A = 3.0, 
< k >= 6 and the network size is 1000. The dynamical 
parameters are a — 1.9, e = 0.9. The length of data used 



- L a =10 4 , e=0.9 

- L d =10 s , c=0.9 
L d =10 s , b=0.1 









FIG. 3: (Color online) Components of principal eigenvector of 
correlation strength matrices of coupled maps on networks vs. 
node degree for ER random graphs, of different data length 
and different coupling strength. The network size is 1000, 
< k >= 6, and a = 1.9. The data has been averaged on 100 
ensembles. 

The Bak-Tang-Wiesenfeld (BTW) sandpile model 
on SF networks has been investigated by K. S. Goh et 
al. 0] , which set the threshold of being toppled of each 
node i be its degree fcj , unlike that of the uniform thresh- 
old height in lattice cases; and proposed a losing proba- 
bility / when the grains are transferring from a toppled 
node to its neighbors, as the sinks or open boundaries in 
lattice cases. The rule can be adopted on arbitrary net- 
works. The correlation in this case is between the top- 
pling events sf in avalanches, sf = 1 if node i toppled in 
the nth avalanche, and if not. It is expected that the 
cooccurrence between nodes in an avalanche could unveil 
the linkages between them. 

The relationship of the components of principal eigen- 
vector of the correlation strength matrix of the cooccur- 
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FIG. 4: (Color online) Components of principal eigenvector of 
correlation strength matrices of the cooccurrence in avalanche 
events on networks vs. node degree for various network mod- 
els. The same panel indicates the same network model as that 
in Fig. 2. 



FIG. 5: (Color online) Components of the principal eigenvec- 
tor of correlation strength matrices of the currents for IFN 
model on networks vs. node degree for various network mod- 
els. The same panel indicates the same network model as that 
in Fig. 2. 



rence in avalanche events on networks and the nodes' de- 
grees is shown in Fig. The network parameters, such 
as the exponent of SF networks, average degree, and net- 
work size, are the same with that in the coupled chaotic 
map cases. The loosing probability is / = 0.001. Again 
the length of the data is = 10 5 , and the correlation be- 
gins after 10000 avalanche events, which is long enough 
to avoid the transitions. Each data has been averaged 
over 300 network ensembles. Also, the positive relation 
between component and degree is clear. 

Coupled integrate-and-fire neuron (IFN) model. The 
dynamics of each node is described by [Tfij ] 



duj(t) 

dt 
dljjt) 

dt 



= -Ui(t)+I +Ii(t), 
T &i j 



where Iq is an external, time-independent and universal 
driven current, e is the coupling strength, and is positive 
in our simulation, A is the adjacent matrix of the back- 
ground network, as defined above, and V is the spiking 
time of node j. When the potential Ui(t) arrives at a 
threshold Ui{t) — 1, it fires and resets to 0. The IFNs 
can be synchronized under certain conditions, which tells 
nothing about the network topology. So we will focus on 
the unsynchronized states, and discuss the current cor- 
relation, that is, the output signal of node i is just its 
current 7j(t). The same relation is shown in Fig. |SJ The 
parameters are: Jo = 1-3, r = 0.5, e = 0.7. The correla- 
tion begins after a transient time of T = 100 time units, 
and performs for another 100 time units. Beside the dis- 
persion of some data, the positive relation is apparent. 

We suggest a method to locate the important or influ- 
ential nodes (namely, the nodes that have large degrees) 
in dynamical systems on complex networks, according to 



the positive relation between the components of the prin- 
cipal eigenvector of the correlation strength matrix and 
the nodes' degrees. First, one can measure the cross- 
correlation matrix from the output signal of each node, 
and the component values of the principal eigenvector of 
the correlation strength matrix could be calculated. The 
degree of a node is represented by its corresponding com- 
ponent value. Thus, the important or influential nodes 
are located as the nodes with larger component values. 
In Fig. HO a typical result of the method is shown for cou- 
pled chaotic maps on a BA network. For other dynamical 
and network models, similar results can be obtained. 

Further more, we define the efficiency of the locat- 
ing method as E(f) — n m (f)/(N ■ f), where N is the 
total number of the nodes, and / is the fraction of se- 
lected nodes with largest degrees, and the n m (f) is num- 
ber of nodes that correctly matched. If the nodes were 
randomly ordered, the fraction that matches with other 
ordering, say, by degree, will be equal to the fraction 
that compared, thus E(f) = f. A typical result by our 
method for IFNs is shown in Fig. [7| The efficiency of 
the method boosts up rapidly as / increases from 0, and 
remains high for almost all the region [0,1]. The high 
efficiency of the method in IFNs model indicates direct 
application potentials, i.e., identifying the most influenc- 
ing neurons in the experiments in neuronal networks |l2| , 
(not necessarily through current correlation). 

In conclusion, we studied coupled chaotic maps, 
avalanche processes, and integrate and fire neuron mod- 
els, covering three types of coupled dynamical systems, 
on various network topologies, whose degree distribution 
varies from Poissonian to scale-free and to exponential, 
and of different local properties. All our results support 
the proposal that the components of the principal eigen- 
vector of the correlation strength matrix have a positive 
relation with the nodes' degrees, thus a representation 
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FIG. 6: (Color online) A direct show of the efficiency of our 
method. Squares: the first 8 nodes with largest degrees; filled 
squares: nodes that are correctly located out by our method; 
empty square: the missed node, which is just the 8'th node 
among the first 8 nodes; small circles: ordinary nodes. 
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FIG. 7: (Color online) Efficiency of the locating method by 
ordering of the component values of the principal eigenvector. 
The parameters are the same as that of Fig. [S] 



of nodes' degrees, the ability to influence others and the 
importance to the system, could be realized by the com- 
ponents, which can be obtained through the dynamical 
processes. Further applications such as immunization of 
internet or contagious disease, identifying pivotal neu- 
rons, etc., can be investigated. 
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